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simulations that confirm the findings of the mean-field theory. 
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I. INTRODUCTION 

Equilibrium and non-equilibrium relaxational dynamics play an important role in many critical phenomena pj. 
Typically, relaxational flows drive the evolution of the system to the equilibria determined by a Lyapunov energy 
functional T that depends on local potential functionals and on the interactions in the system @, Q • The celebrated 
(j) 4 functional is a paradigmatic example of such a potential functional and gives rise to well-known equilibrium 
models such as the so-called model A (a coarse-grained version of the Ising model) , and model B (an archetype of 
phase separation dynamics), among others 0,0)- 

Typical relaxational models describe the flow of a field ifi (t) defined on a d-dimensional square lattice via a Langevin 
equation of the form 

^(*) = -r^^ + r 1 ^(t). (i) 

Here i labels a lattice site, T is a positive constant, ({<^}) = (<Pi, ■ ■ ■ , <Pn) denotes the entire set of fields, and & are 
Gaussian white noises with zero mean and correlation functions 

&(*)&(0> = * 2 M(*-0- ( 2 ) 

Typically, the functional T consists of a local potential V(tp) and an interaction term, 

where K is the coupling coefficient. The left-most sum in Eq. J3J runs over all lattice sites and the right-most sum 
over the 2c? nearest neighbors of a given site i. A keystone in this formalism is the link (the "fluctuation-dissipation 
relation") between the intensity Ta 2 of the fluctuating contribution and the relaxation parameter T through the 
constant a 2 that in equilibrium systems is proportional to the temperature. The relaxation coefficient T affects how 
fast the system relaxes to the global steady state. 

Recent studies have revealed the importance of field-dependent relaxation coefficients T(ip) in the dynamics of these 
systems 0. A generic description for such systems is provided by the Langevin equation 

= -r + [r (^(i))]* (4) 

where T(<p) and [r((/?)] 1 / 2 are both positive. For example, it has been shown that relaxational flows driven by field- 
dependent coefficients may present inverted phase diagrams where ordering effects increase with the intensity of 
the fluctuations. This behavior has been observed in polymer mixtures where spinodal decomposition, i.e., phase 
separation, increases with increasing temperature Q. 

The importance of these flows may be even more pronounced in nonequilibrium systems, and goes well beyond 
the scenarios that lead to inverted phase diagrams. One example is that of pure noise- induced phase transitions Q. 
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Such phase transitions exhibit the striking feature that noise is the crucial element responsible for the appearance of 
ordered phases that disappear in the absence of noise. The first mechanism identified in the literature leading to such 
behavior Q relied on a combination of the so-called Stratonovich drift that arises under this particular interpretation 
of the noise, and collective effects. The Stratonovich drift in these systems leads to opposite dynamical responses 
at short and long time scales. Collective effects generated by the coupling among the field elements can amplify 
short time instabilities (that would die away in the absence of coupling), thus leading to the observed noise-induced 
phase transitions. As a consequence, there was originally a widespread belief that noise-induced phase transitions 
could only be found in systems where there are no noise- induced transitions Q, since the latter transitions occur 
in zero-dimensional (uncoupled) systems. Recent studies involving relaxational flows with field-dependent relaxation 
coefficients have shown otherwise by presenting a system where both a transition (zero-dimensional) and a phase- 
transition (coupled systems) are induced by the same source of noise Q. In fact, the mechanism is not attributable 
to the Stratonovich drift because these transitions occur independently of the noise interpretation [8|. Moreover, the 
same mechanism has been extended to pattern formation phenomena generalizing a previous mechanism based 
on the Stratonovich drift [Icj . 

A thorough understanding of relaxational models driven by field-dependent coefficients is therefore important for 
a number of reasons. They play a relevant role in critical phenomena, they may explain situations where inverted 
phase diagrams are obtained, and they constitute a generalization of the seminal work of Horsthemke and Lefever on 
noise-induced transitions to noise-induced phase transitions in coupled systems. 

Whereas specific non-equilibrium relaxational models driven by field-dependent coefficients have been considered in 
the literature, herein we present a more general analysis of such models, whereby the specific cases considered earlier 
become part of a broad panorama. We ask two questions: (1) What are the circumstances (features of the model, 
values of the control parameters) that lead to purely noise-induced phase transitions? (2) What is the nature of the 
phases that can occur in these systems, and what are the features of the model that determine these phases? In 
answering these questions, we discuss the possible phase diagrams that can be obtained and show that their overall 
structure depends on geometrical properties such as the balance of convexities of the local potentials and of the field- 
dependent coefficients. Furthermore, we show how multistability can be induced by noise. Our point of departure for 
this analysis is the equation obtained by implementing the functional derivative of T indicated in Eq. Q , 

<pm = r fa) (-^^ + + F (fit 2 £<(*), (5) 

C being the discrete version of the diffusion Laplacian operator, 

(ij) 

The paper is organized as follows. In Sec. [H] we construct the mean field approximation to Eq. (J^J and establish 
the different phases that may appear in the model. A discussion of the possible phase boundaries between these 
phases, and the order of the transitions, are presented in Sec. IIIII The structure of the resulting phase diagrams is 
presented in Sec. lIVI as is a specific illustration that corroborates our more general analysis. Sections IIIII and llVl are 
supplemented by an appendix where we show that a particular type of transition, while it may occur, is necessarily 
an isolated point in the phase diagram. In Sec.[V] simulations of the full model corroborate some of our most striking 
results. We conclude in Sec. IVII with a summary and some directions for future research. 



II. PHASE TRANSITIONS: MEAN-FIELD ANALYSIS 



We focus our analysis on systems that may undergo Ising-like phase transitions. A convenient order parameter to 
characterize the phase transitions is akin to the magnetization, 

m=\{<p)\, (7) 

where the brackets indicate both a spatial and a temporal average of the field in the steady state. Ordered states are 
associated with m ^ 0. 

In order that a particular system described by Eq. JSJ experience a phase transition driven by a spontaneous 
symmetry breaking of the order parameter from m — to m ^ 0, the symmetry that leads to m — must be 
embedded in the model. Note that if Eq. J^J) is invariant under the combined transformation 



If <H> -if, 



(8) 
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then m = 0, the symmetric state, is indeed always a solution for the order parameter. Equation (JSJ satisfies the 
required symmetry if V(ip) and T((p) are even functions. We will thus consider this case throughout the paper. 
Furthermore, we also assume with no loss of generality that 

V(0) = 0, 

r(o) = i. (9) 

The exact stationary probability density of Eq. JSJ can be calculated for any noise interpretation, including the Ito 
and the Stratonovich interpretations |8(- However, any further analytic insights require further approximation. We 
implement a mean-field approximation in Eq. by replacing the average value of the fields of the 2d nearest-neighbors 
of any site i by the mean field value (ip), that is, 

(10) 

(ij) 

This procedure, which is equivalent to assuming global coupling rather than nearest neighbor coupling, disregards 
fluctuations of the neighboring sites around the mean value. Since all sites are then equivalent, the lattice index can 
be dropped and the set of field equations reduces to a single equation. However, the unknown mean value of the field 
appears in this equation and must be chosen self -consistently. Thus, we obtain a closed approximate version of the 
problem as expressed in the two equations 



dtp 

(<P) = (<P)p- (12) 
Here (•) stands for a statistical average with respect to the stationary probability density associated with Eq. ((lip. 

p st (<p; (ip)) = N((<p))T(cp)^e-M v (<p)+^(<p)-9r)^ (13) 

N((ip)) is the normalization constant, and a = (a = 1/2) for the Ito (Stratonovich) interpretation of the noise. This 
mean field formulation can not be solved in full generality either, but it does allow some analytic characterization of 
the problem. It is this characterization that we pursue as far as possible. 

Note that the disordered solution (symmetric state) (ip) = always solves Eq. (|12[) . Yet, other solutions such that 
(ip) ^ are also possible. We refer to the latter as ordered solutions. Note that as a consequence of the parity of 
the functions V((p) and T(ip), (ip) p is an odd function of ((p) and therefore, if (ip) is a solution of Eq. (|12fl . then so is 
— (ip). However, both lead to the same value of the order parameter m. 

At this point, we must make a distinction between solutions and phases as determined by the stability of the former. 
We call a disordered (D) phase a macroscopic state where (ip) — is the only stable solution. If the solution (ip) = 
is unstable, and only a solution with (ip) ^ is stable, the phase will be called ordered (O). If (<p) = coexists with 
other stable, but ordered, solutions, the phase will be denoted as multistable (M). 

Since the solution of the self-consistency equation H12|) determines the acceptable values of (cp), it is crucial to 
understand the behavior of (<p) p as a function of (p). Noting that 

— = — ( {( p) p - {v)p ) , (14) 

and applying the generalized Schwarz inequality 

{f{v)) P {g 2 {v)) P >\{f{v)g{v)) P \ 2 (is) 

with f((p) = (p and g(ip) — 1, one concludes that the right hand side of Eq. (|14fl is positive, and therefore (p) p is a 
monotonically increasing function of (ip). Moreover, taking the limit of 

/oc 
dpp n Pst (p;(p)) (16) 
-CO 

as {(p) — > =boo immediately leads to 

lim (p n ) p - (p) n p , (17) 

{tp}— »±oo 
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and, consequently, 



Therefore it follows that 



lim » = 0. (18) 



lim (tp) p <(<p), lim (<p) p > (if), (19) 

{tp) — >oo (tp) — > — oo 

that is, {<p) P necessarily lies below (above) (ip) as (p) goes to plus (minus) infinity. 

Figure illustrates the resulting possible different phases in terms of the possible solutions of the self-consistency 
equation i|12|) . The dashed lines represent (tp), the solid curves (y>) p , and the self-consistency solutions are their points 
of intersection. Our subsequent discussions presume that the system definitely goes to a stable state. We exclude 
"runaway" systems that do not fall into this category. 



III. PHASE BOUNDARIES: SECOND AND FIRST ORDER PHASE TRANSITIONS 

With Fig.^in mind, consider now the possible resulting behaviors of the order parameter as we transition from one 
phase to another by changing a control parameter. Note that the general model 10 (as well as the mean field version 
of the model) depends on only two parameters, the coupling coefficient K and the noise intensity er 2 , so the control 
parameter that characterizes a change from one phase to another could be either of these two (or some combination 
of them). The question then is how the points of intersection in Fig.Hmove as one varies a control parameter that 
takes the system from the behavior shown in one panel to that shown in another. 

Consider the transition as a system moves from the behavior shown in panel (a) of Fig. ^ to that of panel (b) . 
The way this is expected to occur is that the solid curve rotates upward so that in addition to the (ip) = solution, 
another solution emerges at the origin. This second solution then moves upward along the diagonal as the control 
parameter increases. This (a)— >(b) transition is illustrated in the top panel of Fig. [21 where the solid lines represent 
stable solutions and the dotted lines the unstable solutions. The transition between the disordered and ordered phases 
is continuous in the order parameter (second order phase transition) . A transition from panel (b) in Fig. ^ to panel 
(c) would involve the evolution of a kink in the curve that first cuts the diagonal at the origin and then moves upward 
toward the existing cut. Associated with this there is a change in the curvature of (<p) p near the origin. This transition 
is sketched in the middle panel of Fig. [2 where we show the two nonzero solutions moving closer together as the 
control parameter increases. The transition between the ordered and multistable phases is discontinuous (first order 
phase transition) and is expected to exhibit hysteresis. The discontinuity is clearly apparent, for example, in the jump 
from the disordered branch of the multistable phase as one decreases the control parameter and this branch becomes 
unstable. Finally, a transition from panel (a) to panel (c), the bottom panel of Fig. [21 again involves a change in the 
curvature near the origin and the evolution of a kink that first touches the diagonal at a single nonzero value that 
then separates into two as the control parameter increases. The transition between disordered and multistable phases 
is again discontinuous. Note that all of these transitions can also proceed in the opposite direction than illustrated 
here, e.g. a transition from order to disorder would occur as in the top panel but from right to left. 

We stress that the transitions between disordered and ordered phases are continuous in the order parameter (second 
order phase transitions), whereas the transitions from or to multistable phases are discontinuous (first order) and are 
therefore expected to exhibit hysteresis. There may be singular isolated exceptions to this latter conclusion, as 
discussed in the Appendix. 

A more quantitative characterization of the phase transitions is possible for those transitions that involve a change 
in the stability properties of the (p) p = solution, that is, for phase transitions between disordered and ordered 
phases, and between ordered and multistable phases (i.e., the top and middle panels in Fig- EJ - This characterization 
involves the first two nonzero derivatives of (<p) p with respect to (ip) in the vicinity of the origin. The first derivative 
provides information about the slope of {p>) p and the third about the concavity/convexity ({tp) p is an odd function of 
(p) , so even derivatives around the origin vanish) . Transitions between disordered and ordered phases (second order) 
occur when 



d(<p) p 



d(p) 



d{Lpf 



= 1, 

< 0. (20) 

( v )=o 
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FIG. 1: Schematic of the possible solutions for the self-consistency equation 1121 . If the only solution is (<p) — 0, panel (a), the 
system is in a disordered phase. In panel (b), there is another solution in addition to the one at (ip) — 0. Only the nonzero 
solution is stable, so that this represents an ordered phase. In panel (c) there is an even number of nonzero solutions (here, 
two) in addition to the (<p) = solution. The (cp) — solution is stable, the next intersection is unstable, and the third is 
again a stable solution. The system thus exhibits multistability in this case, and involves the coexistence of a disordered and 
an ordered solution. 



Transitions between ordered and multistable phases (first order) occur when 



(v)=o 



> 0. 



(21) 



In turn, these derivatives can be expressed in terms of the cumulants of the probability density pstifl 0). We proceed 
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FIG. 2: Behavior of the order parameter, m, as a function of a control parameter, CP, in the vicinity of a transition point. 
Transitions between disordered and ordered phases are always second order (top panel). As shown in the middle and bottom 
panels, phase transitions from, or to, multistable phases are (with the exception of singular points discussed in the Appendix) 
first order. 



to establish this relation. 

Since (<p) p is an odd function of (ip), its Taylor expansion around the solution (ip) = reads 



(p) p = J2a 2n+1 (p} 2n +\ (22) 



n=0 

where 



00 2 "+V gt (p; (cp)) 
a2n+1 (2n+l)\J_ 00 9(^)2«+i 

Furthermore, a straightforward but tedious calculation leads to the result 



pdtp. (23) 



22n+l f K\ 2U+1 

a-iri+i = -jz — — prr — o C 2n +2-, (24) 



(2n+l)! \a 2 

where C 2n +2 is the (2n + 2) th cumulant of the probability distribution p s t(f, 0). The relation of the cumulants to the 
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statistical moments of the probability distribution is given by 



G. 



2n+2 — 






1 



















1 














<^ 2 >o 
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<^ 6 }o 





O 4 >o 











(25) 



2n+2 



where H 2ll +2 indicates the determinant of the (2n + 2) x (2n + 2) matrix, Q the binomial coefficients, and (-)o stands 
for statistical averages over the probability density p st (y;0). Therefore, around ((f) = the self-consistency equation 
reads, 



00 2 2n + 1 

(</>) = (<p) P = (<p) 



(2n+ 1)! \a 



K 



2n+l 



s c 2n+2 (<py 



and it then follows that 



d(<p) 



2K 



d(<p)* 



-Co 



K 



(26) 

(27) 
(28) 



In summary, the boundaries between disordered and ordered phases and between ordered and multistable phases 
are characterized by the first two nonzero cumulants of the probability distribution p st (tp]0) as follows: 



C 4 < 



c 2 

c 4 > o 



2K 



Second order phase transition (order-disorder boundary) 
First order phase transition (order-multistable boundary). 



(29) 
(30) 



Transitions between disordered and multistable phases can not be characterized in this fashion since they require 
information about (tp) p away from (tp) = 0. 



IV. OVERALL STRUCTURE OF PHASE DIAGRAMS 



Having discussed the possible phase transitions that might be observed in the mean field system l|ll jl -l|12 |l . we now 
ask which particular phases might be present for particular values of the control parameters K and a 2 . In this section 
we present an analytic deduction of the phases present for small and for large values of the coupling coefficient K. 
The behavior for intermediate values must be deduced on the basis of plausibility arguments that we introduce later. 

We start by defining 

IK 

Q K = a x -\ = ^C 2 - 1. (31) 
a" 

This quantity measures the differences in the slopes of (ip) and {ip) p as a function of (<p) near the origin. According 
to the analysis presented in the previous section, Qk > for an ordered phase and Ok < in either disordered or 
multistable regions of the phase diagram. Moreover, the phase boundaries to or from ordered states are given by the 
zeros of 9 ^ . 

Note that Go = —1, that is, (tp) p = in the absence of coupling. Thus, (ip) = is the only possible solution to the 
self-consistency equation (|12|) in the small coupling limit, and the system is disordered in this limit. At sufficiently 
weak coupling the system is therefore always disordered. 
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One can easily check that 



K 



dK 



> 0. 



(32) 



K=0 



Thus, as K grows from zero the system advances toward the ordered phase. This statement does not mean that the 
system will actually enter into the ordered phase as the coupling increases; it simply states the ordering role of weak 
but increasing coupling. 

On the other hand and more interestingly, it is possible to investigate the strong coupling limit as follows. We first 
introduce the convenient notation 



/OO 
V 2n ^)e-^ 2 d^ 
-OO 



(33) 
(34) 



where n ^ 0. We can then write all the non-zero moments of pst{<p', 0) as ((p 2n )o = Iznjlo- In particular, the cumulants 
of interest here can be written as 



h 

V 



h „I 2 



2 - — , U - — - 



Moreover, notice that all the moments can be reduced to the calculation of Iq since 

hn(K,a 2 ) = (-l) n a 2n ' 



,d n I (K,a 2 ) 
dK n 



(35) 



(36) 



A series expansion of Iq useful for large values of K follows from an expansion of $(</s) around (p = 0, which allows us 
to carry out the integral: 



Io(K,a 2 



00 f 



00 m=0 



« $(2m) 7r l/2 

d¥ — m \ 22m 

m=0 



K 



where, 



d z <I>{(p) 



(37) 



(38) 



ip=0 



Introducing Eq. I|37|) into Eq. (|36|l . we obtain the series 

[ m + l/ 2 ] n 



l! 2 2 " 



where 



i-i 



Wn=II + Wo = l" 



(39) 



(40) 



2=0 



The moments of p s t ((/?; 0) then read 
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A more convenient expression for Eq. I|41|l is obtained by performing its Taylor expansion around (a 2 /K) — > 0, 



[l/2]„ + ( ^ 



$(2) 



($(4)$(o) q 5/2 



32 ($(°)) 
(V 2 >) 2 ([3/2]„-[l/2]„) 



([3/2]„ - [1/2]J 
[1/2]J 



(42) 



The first term in the series (Tf2*)l . i.e., up to order (a 2 /K) n , leads to the familiar result of applying the steepest descent 
method 01 to Iq, 



(2 \ n 
^) [1/2L 



(43) 



However, this result is not sufficiently accurate to capture enough of the large- coupling behavior of (<p 2n )o and shed 
light on the behavior of the phase boundaries in that limit. Keeping up to the next order, that is, 



a 

fO»i V K 



[1/2], 



$(2) 

K I 8$(°) 



one finds for the function @k 



_ 1 fa 2 \ 
ex»i- 2 1^)^(0) 



([3/2]„ - [1/2] 



(44) 



(45) 



Note that $ ( -°- ) = 1 [c.f. Eq. ©], and therefore the sign of Qk for large values of the coupling is determined by the 
sign of $( 2 ), 



(2) 



( a _ i)r« - 4^ (2 



(46) 



If <I>( 2 ) > then 6^»i > (ordered phase). On the other hand, if $( 2 ) < then &k^i < (disordered or multistable 
phase). That is, whether or not the system is in an ordered phase depends only on the balance of convexities of the 
local potential and the field-dependent coefficient at the origin. 

Furthermore, since Oif>i = O ((cr 2 /X)) , which vanishes as K increases, one knows that for large coupling the 
system is "near" a phase boundary of an ordered phase. One can gain some insight into the type of transition that 
might be involved by studying the fourth cumulant, C4 [cf. Eqs.J2HJ) and l|3T?l) ]. Using Eq. in Eq. (|33|l and 
recalling that < i>^ - ) = 1, we obtain, 



Ci (K » i) = 1 (!L) _ 3 ($( 2 



(47) 



where, 



$(4) = _if (a _ 1)F (2) r (2) + ^(y{2)) +(a-l)rW 

+3(a-l)(a-2) (r^V - AyW 



(48) 

We thus confirm that, independently of the behavior of the system at intermediate values of the coupling, for large 
coupling the appearance or disappearance of ordered phases as reflected in the sign of the fourth cumulant depends 
on the geometrical properties of V(<p) and of T(tp) around the origin. 

There are basically only two distinct generic types of behavior of these functions around the origin, and therefore 
only four possible combinations. The possible types of functions are shown in Fig. [31 where we have plotted the simple 
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FIG. 3: Generic local potentials Vi(ip) and field-dependent coefficients Ti(tp) as a function of the field ip. The solid lines are for 
i — 1 and the dotted lines for i = 2. The behavior at the origin of the derivatives of V(<p) and T(<p) determine the phase at 
large values of the coupling (see text). 



representative cases |l3| 



and 



2 ' 
4 



2 



(49) 



1 + ip 2 
l + ip 4 ' 
1 



1 + <p 



2 ' 



(50) 



Note that the field-dependent coefficient Ti favors fluctuations around ip = ±1 while ]?2 leads to the largest fluctuations 
around <p = |4|. It is quite straightforward to determine the sign of C4 on the basis of the derivatives of these functions 
(one can use the generic forms l|49|l and (|50|l as a guide) and compile the following table. We emphasize that these 
are results in the strong coupling limit. The entries in the table indicate order (O), disorder {D), and multistability 
(M): 











M 


if cr 2 < cr" 
M if <7 2 > tJ 2 




0\ia A > a' 2 c 
D if cr 2 < cr 2 






Here cr 2 is a critical value of the noise intensity that separates different phases; for the generic models displayed above, 
cr 2 = 1/ (1 — a). Note that the noise interpretation through the value of a simply shifts the critical value of cr 2 . The 



11 



case Vi, T2 with the functions given above has recently been studied in the context of phase transitions and pattern 
formation Hill. Also, in agreement with our general analysis it was noted recently @| that in that particular 
case the noise-induced phase transition is not attributable to the so-called Stratonovich drift, as is the case in other 
noise-induced phenomena [||. 

We have thus arrived at the generic phase structure for the mean field problem (|lll) - (|12|) in the weak coupling limit 
(the system is disordered) and in the strong coupling limit (as shown in the table). For intermediate coupling we are 
not able to provide a general quantitative analysis except to note that if a 1 — > 0, the system is always in a disordered 
(D) phase since the fluctuations are needed to provide the energy to induce symmetry breakings. 

Despite this difficulty, one can introduce compelling arguments to connect the phase behavior that we have estab- 
lished in these limits, and to arrive at a set of full phase diagrams. 

There are several unknown regions connecting various known phases at this point. In the Vi,!?! case we need 
to connect the weak-coupling D phase to the strong-coupling M phase. For the V^Ti combination we require a 
connection between the D phase and an phase if the noise is weak, or to an M phase if it is strong. With the 
Vi , T2 combination the connection needs to be made between the disordered weak-coupling phase and an ordered 
or a disordered strong-coupling phase depending on the noise intensity. And in the Vi , T2 case a connection needs 
to be established from the disordered to the ordered phase. The simplest possible scenarios for connections are the 
following. The simplest connection between disordered phases is simply a disordered phase, i.e., a situation where 
no phase transition occurs at all. A connection between disordered and ordered phases is most straightforwardly 
accomplished through a single second-order phase transition. Finally, for the connection between disordered and 
multistable phases, two different scenarios are most feasible. One possibility is that the connection is mediated 
through an ordered phase, as follows. As noted in Fig. [2 when a multistable region appears from a disordered phase, 
the unstable solution tends at first to move downward as the control parameter increases. If the unstable solution 
eventually vanishes, the disordered phase necessarily becomes unstable and one necessarily enters an ordered phase. 
Such destabilization does not occur if the transition is mediated by and ordered phase. On the contrary, multistable 
phases arising from OM transitions grow more stable as the control parameter increases (see Fig. EJ. In this case, a 
feasible sequence would be of the form DOM. On the other hand, a direct DM transition may also occur, but only if 
the initial vanishing tendency of the unstable solution is stabilized as the coupling increases. 

We can corroborate this scenario by calculating the phase diagrams that are obtained from the mean field approxi- 
mation for the particular functions (|49|l and (|50|l . We numerically solve the self-consistency equation (I12[) and compute 
the boundaries separating different phases. Recall that the noise interpretation simply shifts the transitions points 
but does not change the phase diagram structure. We present the results for a = 0, that is, the Ito interpretation, 
for which a 2 = 1. The results are shown in Fig. 0] 

Note that the actual structure of the phase diagrams at small and large values of the couplings is perfectly captured 
by our analysis. Moreover, the value of the critical noise intensity is a 2 — 1, as predicted. As for the unknown domains 
shown for intermediate coupling, our arguments about the simplest scenarios agree with the mean-field results. For 
example, for the V^Ti combination the appearance of the multistable phase with increasing K is mediated through 
an ordered phase. We point out that for the case V\, Y\, there is a triple point where all phases merge. At this critical 
point a continuous phase transition between disordered and multistable phases occurs. As noted earlier, this behavior 
is singular and isolated. Note also that multistability appears in this case by means of a DOM sequence with increasing 
coupling K for noise intensities to the right of the triple point. However, to the left of the triple point the multistable 
phase arises from a DM sequence. Moreover, above the triple point and with increasing er 2 , multistable phases are 
always destabilized and followed by an ordered phase, as mentioned above. There is evidently an asymptote at a 2 ~ 2 
for the phase boundary separating the disordered and multistable phases. However, this critical noise intensity is not 
captured by our theory since it does not involve slopes and convexities near the origin. Instead, we show below how 
this second critical value of the noise intensity can be calculated by analyzing the zero-dimensional version of the 
problem. 

A number of other striking features of the phase diagrams are noteworthy. For both V\ cases, independently of 
the value of the coupling, the system becomes more ordered as the noise intensity increases. In the case of Ti this 
behavior is associated with the destabilization of multistable phases and in general suggests that the phenomenon of 
the so-called inverted phase diagrams depends mainly on the convexity of the local potential around the origin. It 
is also worth noting the phenomenon of reentrant noise- induced multistability for the case V2,T%. If the coupling is 
greater than K ~ 10, increasing the noise intensity causes a transition from an ordered phase to a multistable phase. 
However a further increase in the noise intensity eventually leads the system back to the ordered phase. This behavior 
resembles the phe nomenon of reentrance as a function of the intensity of the fluctuations in other noise-induced 
phenomena |rj Il0|. However, in the latter the phase changes are from disordered to ordered and, for sufficiently 
intense noise, back to disorder. 

Zero-Dimensional Analysis. Previous studies for the particular case V\, Y2 [H@ have revealed that in the case of 
relaxational-flows with field dependent relaxation coefficients the mechanism responsible for the phase transition is 




FIG. 4: Mean-field phase diagrams as a function of the local potentials and field-dependent kinetic coefficients given by 
Eqs. 14911 and (I5UII . The small open circle in the phase diagram for Ti and Vi where the three phases merge indicates an 
isolated singular critical point (triple point) where a continuous phase transition between disordered and ordered phases occurs 
(see the appendix) . The overall structure of the phase diagrams is in perfect agreement with the schematic structures discussed 
in the text. 



similar to that which drives the noise- induced transition a la Horsthemke-Lefever in zero-dimensional systems Q • In 
the zero-dimensional case, noise-induced transitions are associated with changes in the extrema of the local potential. 
We present an analysis of the zero-dimensional system to compare with some of the results of our mean- field analysis. 
The zero-dimensional version of Eq. jSj reads 

m = re?) (-^p) + n<p)]*m m 

that is, the uncoupled version of our original model JSJ. The stationary probability density now is 

P M = Ne~^ v ^\ (52) 
where N is the normalization constant and V e g(ip) is the effective potential, 

V eS (<P) = V( V ) + g2(1 ~ a) ml». (53) 
The equilibria of the effective potential are given by the condition V^ s ((p*) — 0, that is, 

¥)yV)+ j2(1 2 " a) rV) = o. (54) 

The stability of the equilibria depends on the sign of the second derivative of the potential at the equilibrium points. 
A noise-induced transition occurs when there is a change in the stability of the solution tp*. Therefore, the boundary 
of stability is given by V^((p*) — 0, that is, 

rV)F'V) + a2{1 ~ a) (iv)r'V) - (rv)) 2 ) = o. (55) 
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Note that ip* = is always a solution of Eq. (|54|l and therefore the disordered solution is stable if 



v"(p)+ v V (o) = o. 




(56) 



This equation corresponds exactly to the stability boundary associated with Eq. (|46|) . In other words, the strongly 
coupled system behaves exactly as the uncoupled system. In the cases Vvitp), and Vi(tp), ^(y) the critical value 

of the noise intensity that changes the stability of the disordered solution to an ordered one is exactly as calculated 
in the coupled system, a 2 = 1/(1 — a). 

Moreover, our purpose in analyzing the zero-dimensional system is also to understand the phase boundary that 
separates multistable and disordered phases in the case Vi,ri in the spatially extended problem. Recall that in that 
case the stability of the disordered solution does not change, and therefore we are not able to use our mean field 
analysis near the disordered state to compute phase boundaries. However, we can use Eqs. (|54f) and (|55|) to support 
our numerical findings that indicated that there is a critical noise intensity separating those two phases. While solving 
Eq. H54fl for solutions ip* ^ is rather cumbersome and Eq. (|55|l does not have an analytic solution for those values, 
it is trivial to solve the problem numerically. The result agrees perfectly with our previous findings: there is an 
asymptote at a 2 ~ 2 (Ito) that corresponds to the critical value of the noise intensity separating disordered and 
multistable phases. 



To check the predictions of the mean field theory we present numerical simulations of Eq. JSJ in a two dimensional 
square lattice with nearest neighbor interactions, the Ito interpretation of the noise, and periodic boundary conditions. 
We focus on the case V\,Y\, This is the most interesting, previously unexplored, case: a striking effect of the noise, 
a noise- induced multistable phase and an inverted phase diagram, occur in this case. Note that the case V\ , I^ has 
been studied recently Q , the case Va , T2 presents a phase diagram with a phenomenology similar to the well-known 
model A and the case V^Ti presents as its main feature the same striking phenomenology of noise-induced 

multistability as does the case Vi,ri. 

Figure shows the order parameter to as a function of the noise intensity a 2 for a fixed value of the coupling 
constant, K — 10. The system is seen to explore the three possible phases, disordered, multistable, and ordered, as 
the fluctuations become stronger. Moreover, the system presents an inverted phase diagram where order becomes more 
prominent as the noise is increased. As indicated by the discontinuous behavior of the order parameter, the phase 
transitions are first order in all cases. To detect the multistable phase and the associated hysteresis, we integrated 
Eq. JSJ under two different conditions. The initial conditions and the noise realizations are identical in both cases. A 
difficulty in such simulations is that it takes an inordinately long time to reach a steady state (eventually the system 
leaves any steady state if the system is finite, but this time can be made as long as desired by increasing the size 
of the system) . To overcome this difficulty, in one case we added a very small external field that favors the solution 
to = while in the other we added one that favors an ordered solution. As soon as steady states were reached, the 
external fields were turned-off. The insets of Fig. \5\ show, by means of a density plot, the values of the field in the 
multistable phase for the points A and B for which a 2 w 3. The scale of the density plots is also presented and is the 
same for both insets. The insets highlight the striking feature of the noise-induced multistability and show the two 
possible states within the multistable region. 



We have presented a detailed studied of phase transitions in models with field-dependent relaxation coefficients. 
By means of a mean-field approximation in combination with other analytical techniques and plausibility arguments 
supported by particular examples, we have elucidated the phase diagrams that can be found in such models. We 
stress that our methods can easily be applied to a variety of other systems. Moreover, we have demonstrated that 
disorder-multistability continuous phase transitions are singular points in the phase diagram, and that the phases for 
large values of the coupling are determined by geometrical features of the local potential and of the field-dependent 
coefficient in the vicinity of the origin. We have also showed that the behavior of the extended system at large values 
of the coupling coefficient is equivalent to the behavior of the uncoupled zero-dimensional system. Therefore, the 
mechanism responsible for the phase transitions is similar to the noise-induced transitions a la Horsthcmke-Lefever, 
and is not attributable to the Stratonovich drift together with collective effects involved in other noise-induced 
phenomena 0, IT(i| . Finally, we have performed numerical simulations of a particular case to check the results of the 
mean-field approximation. The numerical results are in qualitative agreement with the theoretical predictions and 



V. 



NUMERICAL SIMULATIONS 
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FIG. 5: Order parameter, m, as a function of the noise intensity, a 2 , for the case VI, Ti and a fixed value of the coupling 
K = 10. The noise term is interpreted in the ltd sense. Note that the system explores the three possible phases (disordered, 
multistable, and ordered) as the intensity of the noise is increased. The insets show the steady states of the field for the points 
A and B indicated in the figure (see text). The noise intensity at these points is a 2 ta 3. 



reproduce the main features of the system, most notably, the occurrence of noise-induced multistability, and of an 
inverted phase diagram indicating that stronger noise induces greater order. 

We envision further modifications of these models that would extend the richness of the observed phenomenology. 
In particular, including other degrees of freedom and considering different couplings increases the complexity of the 
multistability phenomena caused by the noise. Further degrees of freedom could, for example, lead to locking of the 
system either in an oscillatory mode (limit-cycle) or in a stationary state (focus) depending on the initial conditions. 
Considering couplings that favor morphological instabilities could lead to pattern formation determined entirely by 
the initial conditions. Work in these directions is in progress. 
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APPENDIX A: SINGULAR TRANSITION 



In Fig. 01 and the discussion surrounding it we noted that with the exception of singular points in the parameter 
space, phase transitions from disorder to multistability are first order (discontinuous). In this appendix we expand 
on this assertion. By "singular" we mean that if there is a point (a 2 ,K*) in the phase diagram where a continuous 
transition between disorder and multistability exists, then no neighboring points in the phase diagram can present 
the same transition. In other words, continuous disorder-multistability phase transitions are isolated critical points in 
the phase diagram, and there is no plausible continuous function K(o~ 2 ) connecting them. 

Figure shows the behavior of the order parameter m as a function of a control parameter in the vicinity of a 
continuous disorder-multistability phase transition. Note that for that behavior to happen, all three roots of the self- 
consistency equation must vanish exactly at the critical value of the control parameter. Moreover, the convexity of 
(<p) p must also change sign at exactly that value. If such a critical point (crj, if*) exists, it must satisfy the conditions 
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FIG. 6: Schematic behavior of the order parameter as a function of a control parameter in the vicinity of a continuous 
disorder-multistability phase transition. This kind of behavior is singular and isolated in the phase diagram. 



These requirements fulfill the Schwarz inequality, Eq. I|15|). and we must therefore conclude that such a critical point 
is possible. 

Now, assume the existence of the critical point (cr^-ftT,). We investigate the requirements for a neighboring point 
((7* + s a 2 , K* + Ek) to also be associated with a disorder-multistability phase transition. A straightforward calculation 
leads to the following condition to be satisfied by e a 2 , 



Again, the only acceptable solution to this equation is Ek — 0. Therefore, if there exists a critical point in the phase 
diagram where a continuous disorder-multistability phase transition occurs, that point is singular in the sense that it 
is isolated. 
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This can in general only be satisfied if E a 2 = 0. As for ek, it must satisfy 
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